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Abstract 

This paper continues a series of investigations on converging representations for the Riemann Zeta function. 
We generalize some identities which involve Riemann's zeta function, and moreover we give new series and 
integrals for the zeta function. The results originate from attempts to extend the zeta function by classical 
means on the complex plane. This is particularly of interest for representations which converge rapidly in a 
given area of the complex plane, or for the purpose to make error bounds available. 

Keywords: Riemann Zeta function, q-series, Euler-MacLaurin summation, Bernoulli Number and 
Polynomial, Riemann hypothesis. 



1. Introduction and Definitions 

1.1. Extensions To The Complex Plane 

The Riemann Zeta function is classically defined as 



oo 



which is a convergent representation for Re s > 1 . An immediate extension to Re s > is given by the 
variant 



which results from separating summands of the form ^ and (2 .j 1)S . This alternating series converges for 

Re s > and reveals that £ has a pole with residue 1 at s = 1, as 2;=i - f = In 2. 
Probably less well-known is the similar representation 



1 y / 1 1 2 \ 

- ' ~ 1 - 2j \(3i - 2) s + (3/ - l) s ~ (30 1 / 

' oo 



1-3 



1 2 1 



(3i - If (3i) s (3i + 1)' 



where the latter converges as well for Re s > - 1 . 

Other interesting representations, as they seem to not appear in the literature and which extend the zeta 
function even further, include 

• (i = i + Zm ((47V - <w + dw) for Re s > - 1 - 

• (l _ V + h ~ fr)^^ ~ 2i=0 ((6,+l)' ~ (6i+2) s + (6i+3) s _ (6i+4) s + (6i+5) s ) 3I1< ^ 

• (4-|-^-^)fW=4- ^-Z ! -=i(^-^ + ^-(^ + (^)forRe,>-3; 
they are all verified in a similar way as the initial formula. 
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1 INTRODUCTION AND DEFINITIONS 



1.2 Euler-Maclaurin 



-dx (2) 



1.2. Euler-Maclaurin 

To further extend the Zeta function to complex arguments with even smaller real part one may apply 
Euler-Maclaurin's summation formula, which states that 

i=\ "' 1 j=i J- J 1 

f is a function with (continuous) derivatives /® up to order A, {x} — x - \_x\ is the fractional part of x, Bi 
the Bernoulli number and B, (x) = Zy=o ij)^i-j x ^ me Bernoulli polynomial (cf. [AS64]). 
The statement may be applied to (1), giving thus after some rearrangements the identity 

((«) = «!r(z-«+D = n ' f^i ' ' ' ^T 1 is the (generalized) binomial coefficient.) 

The error term dx often is inconvenient, as it involves an integral of a function, which is not 

smooth enough. We shall elaborate different expressions below, replacing this error term typically by other 
series or by different integrals; as an example recall the integral 

f W = 1 + — f { k+ S " 2 W + — f — - t -t k - l e-'t^dt (3) 
(cf. [CS10]or[HC01]). 

This construction plays even a prominent role in the famous algorithm of Odlyzko and Schonhage's to 
evaluate ( on the critical strip, so some remarks seem to be appropriate (cf. [OS88], [C092] and [PR07]). 
Remark 1. Identity (2) is frequently given in the form 



fw-z'-z-Mr^-c-Wrl- /™<* 

which is equivalent because - (except for B\ = — 5) and ^T 1 ) = (-l)- 7 ( s+ j 2 )- 

Remark 2. As is a polynomial, B& ({x}) is uniformly bounded for and it even holds that \B& (x)\ < 

(i-2' 2 ^j(2/r) A for x 6 [0 ' 1] (cf [ Leh4 °])- ™ s allows to conclude that 1 

h CD 



i=l /=() N ' 



and whence, f is analytic in the entire complex plane except for s = 1, where is a pole with residue 1. The 
argument, however, can be refined for even more explicit error bounds and rates (cf. [C092]). But as a 
general patter the expression converge, provided that Re s > 1 - A. 

Remark 3. Notice as well that ( A s ) = for s a negative integer and A big enough, which enables to recover 

£ (-71) = — 

n + 1 

forn e {1,2,...}. 

Remark 4. The partial series is frequently denoted 



H n (s):= J i 



which is the generalized harmonic number of order n of s. 
Remark 5. For future reference recall that 



i/„ ( -fr)=x ^+ i( 7j;-^' ( 4) 

1=1 

in the notation just introduced for natural numbers k e {0, 1,2, . . . } - this is sometimes referred to as 
Faulhaber's formula. 

'Note again that By+i = 0, so the order of convergence in the next statement actually is O ( 
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2. Results And Extension To Various Directions 

It turns out that a lot of different relations for the Riemann Zeta function again involve the expression 

('"0 

T,%o ^j~r ■ This * s natural ' as we have seen that this expression describes the asymptotic behavior of 

So does the following extension of an identity which is sometimes referred to as Stark's formula; an 
extension to Hurwitz Zeta function is known under the name Stark-Keiper formula: 

Theorem 6 (Generalization of Stark's formula). For s + 1 and any A e {1,2, . . . }, 

Remark 7. The latter series converges for all s ^ 1. As ( (s + (i > A) are easily available even by means 
of formula (1) provided that Re s > 1 - A, the statement makes available by different choices of A for 
Re i < 1. 

The particular statement for A = 1 (sometimes called Stark's formula, although it seems it was published 
first by Landau) is (for example) contained in [Tit86]. 

Proof. Notice first that 



i-— TT^l 1 "}^-— FAfMt'-i 

s - 1 n s+ '- 1 \ i J s-l^n 5 ' 1 \\ n 

n=2 i=2 x ' n=2 \\ ' 



1-^ 



* - 1 *-i (n - IY « « s - 1 

n—l v 7 



s- r 

which is the statement for A = 1 . 

The general result follows by induction on A: To this end assume that (5) holds true and replace s in (6) 
by s + A. Subtracting appropriate multiples of (6) from (5) to eliminate £(s + A) reveals the assertion for 
A + 1 , where additionally the identity 

A 



k=0 

(A > 1) is involved. □ 



Theorem 8. For A e {1,2,...} and s + 1 the relation 

k=A N ' j=0 • ■ 



holds true. 



Remark 9. The latter series converges for all s + 1. As above, f [s + i) is easily available by means of (1) 
for A > 1 - Re s, the statement whence again makes f (s) accessible for arbitrary Re s < 1. 



Proof. The Statement for A = 1 is 

's + k- l\£(s+k) 



k / 2* 



(8) 



2 RESULTS AND EXTENSION TO VARIOUS DIRECTIONS 



and sometimes called Ramaswami's formula. This is verified straight forward, as 

Y (s + k- 1U0 + &) y 1 ^(s + k-l\(l° 
lj\ k J 2* -LTsIj\ k j\2i 



/=1 N ' 1=1 



(2 s -D a^). 



As for the general statement one may replace s in (8) by s + A and substitute this for the first summand in 
(7). This burns down then to the identity 



A-1 



y Bj y (i) ill = y Bj ( j ) 

Zj 2 s - 2 [ -> Zj 2 J+A - 2 i + 1 Zj 2 s - 2'-> A + 1 ' 

7=0 i=7 7=0 



which holds true because 



y B > 1 y m + 1 \ J_ 

Zj 2 s - 2 l -J 2 !+A - 2 \ //\ ; /i + 1 



7=0 <=7 

A-1 



v g, rr) 2 w- 2 i-A 

Zj 2 s - 2'-> A + 1 2 s - 2 X " A 

_ a y' s .m/_! L_ 

Z-i ; A+1\2 s -2'-> 2 ¥ - 2 : - A 



7=0 

A 



.V Bj ft') 

Zj 2* - 2 l ~i A + 1 

7=0 

andi; A : : 1 B j ( A ; 1 ) = -(A+l)B A . □ 



Remark 10. Equation (8) is similar to other representations found in [Sri03]. 

The ingredients from the previous results can be processed together to give a representation of (2) which 
does not involve integrals: It is the compelling advantage of the following representation over (2) and (3) 
to replace the integral by a usual sum. As the error term can be evaluated precisely, this makes ( available 
without computing the limit and letting n — > oo in the respective sum. 



Theorem 11. For n e {1, 3, ...},/> e {2, 3, ... } and A e {0, 1, 2, ... } the representation 



A-1 ( s+ *~ 2 ) 



1 n 



s+k-l 



J y C*?l) y /' + A B (nip - W +1 ~ J 



(9) 



Zio/cfc frwe, which represents a converging series at least for Re s > 1 - A. 



Remark 12. Notice that (9) even holds for n — 1. 
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Remark 13. The order of summation must not be changed in (9): As formally 

(< +r l) £ + jv (n(p _ 1))i+W ^ 

£u p s_ p i-j£u\ j j (i + l)(np y K 
_ y »J (T) 

which characterizes the asymptotic behavior, but this series does nof converge. 

Remark 14. For p — » 1 and (formally) interchanging this limit with the sum the (non-converging) identity 

1 °° CD B k 



is obtained; and for p — > oo, the same argument gives 

£ CD 



;=1 J A=0 



Proof. As for the proof notice that 



J p' 5 - p x -J 



_ A ^Cr)^_ Ly g/ 1 y( /+r V% LV 

£j s - 1 n^-i B ' 2j (n <p _ !));-! pS - p i-J Zj i + l \ j j \ p ) 

+ 1 y y ( f + % (»(p-i)V +w 

» i ^(/ + i)(«p) , ^\ j J 1 P °-p i -j 

Next we have for > 1 the geometric series 

1 (n(p-l)) M -j (p-l) M - j 1 



_ (p - 1)' +W y 1 y (n(p-l)p k ) 

and thus 

w _ A y'Cr)_^ yji i ^yMfi+i)/, iy 

Zj s _ 1 „«^i n s + j-i {p _ pS _ p i-j La i + 1 \ j J \ pj 

fi+s— 1\ ; . 

i=0 fc=0 ^"f J /=0 v J ' 

The inner series has an explicit evaluation, as 
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For this and the fact that B2J+1 = (except for B\ = — |) the first two sums collapse to \, and thus 
Next observe that the inner sum is just 

7=0 V J ' 7=0 

(cf. (4)), thus 



(^l + y'^-jy 1 . . y / 

1 // + 5 _ A \ n(p-l)p l -\ , 



(«p i+I ) S 

In this situation we may evaluate the power series and rearrange the resulting terms as 



l l n(p-l)p*-l 

( * ) "^ + S(»P t+i ) s ' § 



n(p-l)p*-l 



1 



V 1 n / 



j mp-L)jr-i , | ys 

^ + § § Up* +1 -;7 

k=0 j=l+np t W ' /'=« 7 

B-l , 



7=1 



which finally completes the proof. 

To analyze convergence it should be mentioned that 



/j-l+A 

J_ y v j ' ) V [' + 



(n(p-l)) i+l -j 



1 p s - p l -j 



--y 

n s 

1=A 



1 v / s - 1 + i 



S 



— — J y /-Lyf^l.^.u-D/f 1 ' 1 



i A-l /. , 
1 T-i /I + 1 



Recall from (2) and (4) that Z>! / = O(f^), and further that 



7=1 7=0 V 7 
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Whence convergence can be described by 

j (n(p-Y)p l -l 



=0 



((n(p-l)p k ) 



\i+l-A \ 



=(9 



iV i 



PI p 



k(s-l+&) 



which shows that the series converge at least for Re s > 1 - A. 



□ 



Woon's formula (cf. [BBCOO, Woo98]) is to some (minor) extend close to the representation above, it 
involves a double sum as our representation (9). 



3. Parametrization by Complete Elliptic Integrals 

Riemann's Zeta function has the representation 



^ _1 -(# 3 (0;i0- l)dr, 
o 2 



(10) 



(t) = #3 (0; it) = 2nez <f ' — 1 + 2 Ym=\ <t being Jacobi's elliptic theta function (we choose q = e~"" here 
for convenience; [Edw74] uses ¥ for x ¥(i) = \ (# 3 (0; it) - 1) = £„ =1 e~™)- The latter representation (10) 
deduces immediately from 

rJ dl 



and holds true for Re s > 1 . In order to get a formula for all complex numbers s one may apply integration 
by parts twice for 

rff) r 

sCs-lK(s)-^ = / f5 _1 0(r)dr, 

7!" 2 J0 



(11) 



where 



6 (t) = 2t 2 6" (?) + 3? 6»' (f) = U (n 2 ntf - 6n 2 K?j e^ 2 " 



It can be proved directly that 8 inherits the property 



'(f) 



17=1 



1 - 1 



(12) 



from 0, but additionally 



= lim P k) (x) = limg w (.x) 



for all derivatives including the function 6 itself (k = 0, 1 ... ) holds true, such that the next integral (11) 
defines an analytic function for all values of s e C. 

The symmetric version (for eJ0(e x ) = e~4 0(e~ x ), again by (12)) of (11) reads 



r (0 



;f04) .gf0( e *)d*, 



(13) 



and it is an interesting observation that this is the continuous Fourier transform of x i-> cSCe*): 

r(t) 



s(a-iK(j)-*f 

7T2 



V2^.^( e ^(^))(i(,-i)), 
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where T (/) (k) = -j== e lkx f(x)dx is the Fourier transform of /. 

The property (12) may be used to easily deduce Riemann's functional equation. The function decays 
rather quickly to 0, for e*6(e x ) ~ 4-7i 2 ei x e~ ne 1 as x tends to infinity. 

A good proxy of (11) in terms of a polynomial in s with convergence as n — > oo is given by the 
asymptotic properties of the Hermite polynomial H n (cf. [P6127]) and Laguerre polynomial L„, 



(-D" ■ nl 
(2n)\ 



i \ 



("?) 
1 



4i yfrij 

rH) 



r("D 



X 

n 

£_ 

16/7 



e*6(_e x )dx 



(14) 



4/ 



eW(e x ) Ax 



l \2\ 



e*6(e x ) dx 



r(f) 



The symmetric integrand is maintained by the following form in terms of complete Elliptic Integrals of 
the first and second kind. 



Theorem 15. For any s e 



s(a-l)f(j) 



r(f) 



K(l -m) 
K(m) 



j( s -D 



t/ (m) dm, 



(15) 



where U is the symmetric function 

(K(w)K(l -m))3 



U (m) 



7r V27T m ( 1 — m) 



3(1 -m)-K(m)E(l - m) \ 
+3w ■ E(ra)K(l - m) 
-3 -E(m)E(l -m) 
-mil — m) ■ K (m) K (1 - m) , 



Here, K(m) = L 2 — p 



== = I 2^1 |, l,mj is f/ze complete elliptic integral of the first kind, one/ 

E (m) = j 2 Vl - m sin 2 xdx = | 2^1 (5, -5, 1, m\ the complete elliptic integral of the second kind 2 , which 
describes the circumference of the ellipse with eccentricity yfm. 

We shall use the abbreviations K = K (m), K' = K (1 - m), E = E (m) and E' = E (1 - m) in the sequel. 
A power series analysis of U gives that 



U(m) = 



1 

4^ 











-I) 




\ m 







+ 0(1) 



^Elliptic integrals are often defined with the parameter m = k 1 instead; in the present context, for symmetry, it is more convenient 
to use the parameter m. 



8 



3 PARAMETRIZATION B Y COMPLETE ELLIPTIC INTEGRALS 



and 

K(l-m) 1 16 

v , : = - In— +o(D. 

K (m) 7r m 

As 

In— I dm = 16 ■T\a+ l,ln — ) a x|ln — 

o \ m) \ x j \ x 



1 a 
1 + — - + .. 

ln^ 



is integrable for any a with a result expressed by the upper incomplete Gamma function F. It is thus evident 
that the representation (15) converges for all values of s. 

K(l-m) 

Proof. For the nome q (m) :- e K| »» the q-series identity 

l+2^(mr 2 = ^K(m) (16) 

n=l 

holds true. Observe that 

d , d K(l-m) EK'+KE'-KK' tt 2 

dm 5 ^ dm K(m) 2m(l-m)K 2 4m(l-m)K 2 

due to the well-known derivatives for -^-K = and -^-E = ip^ and Legendre's relation (cf. 

dm 2m(\-m) am 2m ° v 

[EMOF53] or [GROO]). 

Faking the derivative of the formula (16) reveals the closed form 



2 q (m)" 2 n 2 7T = ^K(m)|- (E - (1 - m) K) 



«=1 

after simplification, and proceeding in the same way further 



2 /2 K 2 

q (mf (n 2 Jr) = J -K (m)^ (3E 2 - 6 (1 - m) EK + (3 - m) (1 - m) K 2 ) . 

«=i 

Combining the latter equalities and after further simplifications involving Fegendre's relation again one 
obtains 

t (m)~* ■ 6 (t (m)) ■ t' (m) = U (m) , 
where t (m) := - l - log q(m) = Whence 

r(f) r , 

s(s-l)((s)-^ = / f;-'g(f)df 
tt 2 Jo 

= / t (m) 5" 1 (f (m)) f ' (m) dm = / t(ni)^~ t (m)~ *6{t (m)) f ' (m) dm 
Jo Jo 

7 K(l-m) ^H) 

Tr/ ; (>«) dm, 

o \ K ( w ) / 

which is the desired assertion. □ 



Fhe results found can be combined in the following way. 
Recall the inverse elliptic nome function q~ l (q) := 
this choice x^ 1 (x) - q~ x {e~ ne j is the inverse function of x(m) : = log 



Recall the inverse elliptic nome function q' 1 (q) := (cf. [WW62, Chapter XXI, p. 486] ) 3 ; with 

m)_ 
m) ' 



3 This is called the problem of inversion, cf. [EMOT53, volume II, p. 362]. 
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This can be used to rewrite (13) in different variants as 



*(*-!)?(*)■ 



-dx 



(*4> . # 



x' (jc 1 (x)) 

V # 2 (e _;re ") 



i?3 {e~ m ' x ) 
dt 



dx 



'/ logo\ 2 1 - 
■ ' —U 

qn 



dq 



where 



U(m)=4 



(K(m)K(l -m))J 
n 2 yf2n 



3(1 -w)-K(m)E(l -m) \ 
+3m • E(m)K(l - m) 
-3 -E(m)E(l -m) 
—m (1 - m) ■ K (to) K ( 1 — m) , 



Introducing B (to) := E — ^— = Jq 2 . cos = = =d<ft = | 2-F i (5, \, 2, to) this function rewrites as 

C7 (m) = V8^ W(1 7 OT) (K'K)' (2K'K - 3B'B) . 
Comparing the identity with (13) we finally obtain that t J#(f) = u{i^pJ^ j = j' 1 



ei0(e x ) = t/ 



#3 {e-" e ") 



by symmetry and (12). The inversion moreover has the well-known (cf. [EMOT53, Volume II, p. 362], 

correcting the misprints there) expansion = (y\h=o j+f^+i") as an mnn i te product. 
The power series expansion 



r (f) 

s(*-iK(*)^ 



x 2k ■ e i e{e x )dx 



i 

k=0 



(2k)\A k 



x 2k -U 



&3 {e-" e ") 



*\\ 4 \ 



dx 



for even and positive coefficients then is just a by-product. The convergence of (14), however, is faster. 



4. Summary And Acknowledgment 

4.1. Summary 

Various algorithms to evaluate the Riemann f function start with the Euler-Maclaurin summation for- 
mula, as it properly describes the asymptotics of the series representation. The closed forms available for 
the error bound typically involve integrals. The method proposed here generalizes other methods and gives 
the precise error bound as sum. 

4.2. Acknowledgment 

We would like to thank the reviewers for their significant contribution to improve the paper. 
The plots have been produced by use of Mathematica. 
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